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INTEGRAL monitoring of the Galactic Plane in the last 5 years revealed a new subclass of High 
Mass X-ray Binaries (HMXBs), the Supergiant Fast X-ray Transients (SFXTs). They display 
flares lasting from minutes to hours, with peak luminosity of 10^^ — 10^^^ erg i^' and a frequent 
long term flaring activity reaching an X-ray luminosity of 10^^ — 10^^ erg as recently de- 
tected by the Swift satellite. The quiescent level is around 10^^ erg i^'. We performed a system- 
atic re-analysis of archival INTEGRAL data of four SFXTs: IGR J16479-4514, XTE J1739-302, 
IGR J17544-2619, IGR J18410-0535. This led to the discovery of previously unnoticed outbursts 
from IGR J16479-4514 and IGR J17544-2619. We discuss these results in the framework of the 
different structure of the supergiant wind proposed to explain the outburst from this new class of 
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IGR J16479-4514 


N. 


Start Time 


Peak flux (20 - 


60 keV) 


det. 




(UTC) 


erg cm^^ 




significance 


1 


2003 Feb 02, 19:06 


5.9 X 10- 


-10 


5.95 


2 


2004 Aug 09, 02:24 


1.3 X 10" 


-10 


5.46 


3 


2004 Aug 20, 07:26 


7.0 X 10" 


-10 


5.38 


4 


2004 Aug 20, 12:04 


1.2 X 10 


-9 


7.28 


5 


2004 Sep 10,01:12 


5.8 X 10" 


-10 


6.04 


IGR J17544-2619 


1 


2004 Feb 27, 14:16 


1.0 X 10 


-9 


6.29 


2 


2006 Sep 20, 10:00 


2.3 X 10 


-9 


8.39 



Table 1: Summary of the new flai-es of IGR J16479-4514 and IGR J17544-2619 discovered in this work. 
The peak flux is calculated on a time interval of 200 s. 



1. Introduction 

The INTEGRAL monitoring of the Galactic Plane led to the discovery of many new High Mass 
X-ray Binaries (HMXBs), and in particular of a new class of X-ray transients with OB supergiants, 
the Supergiant Fast X-ray Transients (SFXTs) ([p|]; [|l6|]). SFXTs are characterized by fast X-ray 
flares with peak luminosity of 10^^ — 10^' erg s^^, and a frequent long-term flaring activity with a 



level of X-ray luminosity of 10^^ — 10^^ erg s ^ [|19|]. 



The main hypotheses proposed to explain the SFXTs behaviour are based on the structure of 



the supergiant wind (see Sidoli 2008 for a recent review []20|]): flares could be due to accretion of 



clumps in a clumpy spherical wind (in't Zand 2005 ^) or to enhancement of accretion when the 



neutron star crosses an equatorial wind component (Sidoli et al. 2007 ||18|]). 

Another model involves the presence of magnetars in SFXTs, and the short outbursts result 
from transition across the magnetic and centrifugal barriers (Bozzo et al. 2008 and references 
therein [|]). 

Here we report the results of a systematic analysis of public archival INTEGRAL data of 4 
SFXTs: IGR J16479-4514, XTE J1739-302, IGR J17544-2619, IGR J18410-0535 from 2003 to 
2007. We compare the observational results on these 4 SFXTs with the predictions of the spherical 
clumpy stellar wind model that we have developed. 



2. Observations, data analysis and results 

Using OSA 7.0, we analysed 7328 ScWs (IBIS/ISGRI) ranging from 2003 to 2007, corre- 
sponding to a total exposure time of 15 Ms for the 4 SFXTs. We extracted the lightcurves in the 
energy range 20 - 60 keV and we found 31 flares for IGR J16479-4514, 43 flares for XTE J1739- 
302, 16 flares for IGR J17544-2619 and 3 flares for IGR J18410-0535. Besides the outbursts al- 
ready reported in the literature, we discovered new flares, listed in Table ([l|). For each outburst, we 
have extracted an IBIS/ISGRI spectrum in the energy range 22 — 100 keV and we have performed 
a fit with different spectral models: power law, bremsstrahlung, Comptonization model (COMPTT 



2 



Supergiant Fast X-ray Transients: interpretation of archival INTEGRAL data 



L. Ducci 



(days) (keV) (13 d.o.f.) 

IGRJ16479-4514 47 ThO^^ Tm 

XTEJ1739-302 118 22.2+',';^ 1.71 

IGR J17544-2619 118 ^OA^qI 1.56 

IGRJ18410-0535 35 26.4 +Jj 0.94 



Table 2: Best fit parameters of the average spectra of IGR J 16479-45 14, XTE J 1739-302, IGR J 17544-26 19 
and IGR J18410-0535 with a bremsstrahlung model (IBIS/ISGRI). Texp is the net exposure time. 



in XSPEC). Since there was no evidence of a spectral difference between the flares, within the 
uncertainties, we extracted a total flare spectrum for each of the 4 SFXTs. We obtained the best fit 
with a bremsstrahlung model (Figure [l|), with the parameters reported in Table (^. 



3. Clumpy stellar winds 



We developed a spherical clumpy stellar wind model in order to compare predictions of this 
model with the behaviour of the SFXTs. We assumed that the OB type supergiant is surrounded by 
a clumpy and spherically symmetric wind, where the clump formation rate distribution is: 



N = kM^i^ clump s ^ 



(3.1) 



where is the mass of the clump, and N is the rate of clumps with mass Md emitted by the star. 
We assumed spherical clumps, with radii Rd, then we introduce a power law distribution of radii 

y -I 

Nm °= R[.i clump s (3.2) 



We assumed that the total mass loss rate of the supergiant is given by: 

Mtot=Mcl+M„ind 



(3.3) 



where Md is the component of mass loss rate due to the clumps and M^iW is the mass loss rate in 
the form of the tenuous inter-clump medium. 

From equations (3J) and (33) we obtain that the total mass loss rate of the supergiant, Mtot, 
is given by: 



Mtot = M„ind + / NdM = M,,ind + 



Mb 



M„ 



(3.4) 



In equation ( |3.4| ) Ma and Mb are the limits of the clump mass range. Assuming that Mtot is known, 

fMtot 



from Equation (34) we can determine k: 

fMto, 



(3.5) 



where / = Md/^tot is the fraction of wind mass contained in the clumps. 
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Figure 1: Bremsstrahlung fit and residuals (in units of standard deviations) to the average IBIS/ISGRI 
specti-a of IGR J16479-4514, XTE J1739-302, IGR J17544-2619, IGR J18410-0535. 



Clumps are driven radially outward by transfer of momentum from UV photons to the ions 
of the wind via absorption or scattering in spectral lines From spectroscopic observations of 
O stars, Lepine & Moffat (2008) suggest that clumps have the same velocity law of a smooth 
stellar wind. We can then assume for the clump velocity profile v(r): 



v{r) 



1-0.9983- 



(3.6) 



where Voo is the terminal wind speed, Rqb is the radius of the supergiant, 0.9983 is a dimensionless 
parameter which ensures that v{Rob) ~ 10 km s^^, and j8 = 0.8 is a constant ([|8|]; [^]). 

The clump size is determined by the balance pressure equation. Following Lucy & White 
(1980) 1 10] and Howk et al. (2000) ||5|], we find that the clump size increases at larger distances 
from the supergiant star (see our derivation in Romano et al. (2008a) [pj|]) with the equation: 



Rci{r) =Rci{Rob) 



RqB^oJ 



(3.7) 



where vq = v{Rob) is the initial velocity of the clump at the surface of the supergiant. We found, 
for each mass of the clump, the upper-Umit and the lower-limit for the clump radius by means of 
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two conditions. For any given mass there is a minimum radius below which the clump is optically 
thick in the UV resonance lines. Then the gravity dominates over the radiative force due to the 
line scattering and the clump tends to fall back onto the supergiant. In order to be accreted by the 
compact object of the binary system, the clump must escape from the OB supergiant. Then, from 
the theory of radiatively driven stellar wind of Castor et al. (1975) and the clump model of 
Howk et al. (2000) [f], we found for each mass of the clump the lower-limit for its radius. With 
regard to the second condition, the clump is defined as a density enhancement in the smooth stellar 
wind. Then, for each mass of the clump, there exists an upper-limit for the clump radius: clumps 
with larger radii would be less dense than the smooth stellar wind (inter-clump medium), in contrast 
with the clump definition. The functions of the upper-limit and lower-limit for the clump radius are 
represented in Figure 




2 - 



10 



15 20 
loglO(Mcl) (g) 



25 



Figure 2: Functions of the upper-limit (dotted line) and lower-limit (solid line) of the clump radius at 
f = Rob- In order to obtain these functions we have assumed the following parameters for the supergiant 
wind: Mqb = 30 M©, Rqe = 23.8 Rq, Voo = 1700 km s"', j3 = 0.8, vq = 10 km s-\ C = 2, 7 = -1.5, 
Ma = 5x 10''^ g, Mb = 1022 g, M,,IM,,M = 0.85. 



In conclusion, it is possible to derive the properties of the supergiant stellar wind assuming the 
clump mass distribution (3J), the clump radii distribution (3^), the expansion law of the clump 
( 0^ ) and the upper-limit and lower-limit laws for the clump radius (see Figure 

Assuming the following parameters for the supergiant wind Mqb = 30 Mq, Rqb = 23.8 



,22 



v<» = 1700 km s-\ p = 0.8, vq = 10 km s-\ C = 2, 7 = -1.5, M« = 5 x 10^*^ g Mh = 10 
Mci/M.„ind = 0.85, we calculated different histograms of the flare luminosities for different orbital 
periods and eccentricities of the binary systems (see Figure 



4. Comparison of observed flares with theory 

We have appUed our spherical clumpy wind model to the INTEGRAL archival observations of 
IGR JI6479-45I4 and XTE JI739-302, for which we found a significant number of flares. Our 
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Figure 3: Theoretical histograms of the flare luminosities for different orbital periods and eccentricities of 
the binary systems. The time interval for each histogram corresponds to the orbital period. 



preliminary results are shown in Figure (^. For these two sources we obtained the calculated 
histograms assuming the parameters for the supergiant winds reported in Section (|3|), in a time 
interval equal to the exposure time of IGRJ16479-4514 andXTE J1739-302 (seeTable|). With the 
spectral parameters of Table (@) we obtained the flare luminosities in the energy range 1 — 200 keV, 
for a distance d = 4.9 kpc for IGR J16479-4514, and a distance d = 2.7 kpc for XTE J1739-302 
[12]. The comparison between observed distributions and the calculated histograms reported in 
Figure (|) is limited by the fact that INTEGRAL can detect only the "high luminosity tail" of the 
flare distribution (see Figures ^ and Q). 

The flare luminosity distributions are well reproduced with Port = 2>0d and e = 0.4 for IGR J16479- 
4514, and Port = 10 d and e = 0.4 for XTE J1739-302. With the Kolmogorov-Smirnov test we 
found a probability of 70.6% that the observed and calculated histograms for IGR J16479-4514 
have the same distribution. In the case of XTE J1739-302 this probability is 7.3%. Moreover, 
the observed flare duration is well reproduced by the calculated flare duration (~ 10^ s). The low 
probability obtained with the Kolmogorov-Smirnov test for XTE J1739-302 indicates that for this 
source we need an higher number of flares for the comparison between the calculated histogram 
and the flare luminosity distribution. We note that the estimated orbital period for XTE J1739-302 
is different from that obtained by Blay et al. (2008) [0] (Port ~ 8 d), and that an orbital period 
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of 70 (i with an eccentricity of 0.4 could explain the absence of flares during the whole Key 
Programme 2 period, lasting 32 days, reported by Blay et al. (2008). 




Port = 30d,e = 0.4 
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Figure 4: Comparison of histograms of observed flares with theory for IGR J16479-4514 and XTE J1739- 
302. Dashed lines mark our calculated histograms, solid lines mark observed luminosity distributions. 

In this framework we can study the properties of the supergiant winds of the SFXTs by means 
of observable flare luminosities, flare frequency, flare durations, and quiescent level of the source. 

These are preliminary results on the comparison between the spherical clumpy wind model 
and X-ray observations. Our future work will deal with anisotropic clumpy winds, where there is a 
denser clumpy wind component in the form of an equatorial wind disk around the supergiant star, 
and a polar spherically symmetric clumpy wind component, as proposed by Sidoli et al. (2007) 



[]18|]. We will also compare this model with the data obtained from the still on-going monitoring 
with Swift of these 4 SFXTs [|l|]; [0]. 
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